Excess noise depletion of a Bose-Einstein condensate in an optical cavity 
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Quantum fluctuations of a cavity field coupled into the motion of ultracold bosons can be strongly 
amplified by a mechanism analogous to the Petermann excess noise factor in lasers with unstable 
cavities. For a Bose-Einstein condensate in a stable optical resonator, the excess noise effect amounts 
to a significant depletion on long timescales. 
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A Bose-Einstein condensate (BEC) in a high- finesse 
optical resonator realizes a peculiar hybrid system of 
strongly coupled matter and light waves. The conden- 
sate with its large optical density creates a substantial 
refractive medium for the light field, meanwhile the cav- 
ity field exerts significant mechanical forces on the mo- 
tion of the atom cloud. The study of quantum proper- 
ties arising from this nonlinear dynamics has just begun 
[2 |5J |3J |31 [5] , accompanied by an experimental progress 
in loading ultracold atoms in the tiny volume of a high- 
finesse microcavity [SHU]- 

Many of the properties of a degenerate quantum 
gas can be satisfactorily accounted for by assuming 
a 'condensate wavefunction' which obeys a non-linear 
Schrodinger equation. Although the atom-atom colli- 
sions constantly kick out atoms from the one-particle 
ground state, even at zero temperature, the fraction of 
noncondensed atoms due to s-wave scattering is usually 
negligible. Large depletion of a BEC can indicate some 
inherent many-body interaction effect. In case of colli- 
sions, such observation requires large local densities. An 
example is the strong depletion demonstrated as a pre- 
cursor of the phase transition from superfluid to Mott 
state in an optical lattice [TUITS]. In this Letter we will 
show that a condensate dispersively coupled to the ra- 
diation field in a cavity can also be subjected to strong 
depletion. In this case the many-body effect originates 
from the long-range atom-atom coupling mediated by the 
cavity mode, therefore a significant amount of depletion 
may occur at low density. 

We consider an open system, a pumped and lossy cav- 
ity field interacting with a BEC inside. This system has 
a dynamical equilibrium which self-consistently defines 
the condensate wavefunction. Depletion of the BEC is 
defined as the total occupancy of motional modes other 
than the macroscopically populated BEC wavefunction. 
For the atoms in a cavity, the quantum fluctuations of the 
radiation field is a possible source of depletion. However, 
such a fluctuating dipole potential alone cannot imply the 
unexpectedly large depletion in the dispersive interaction 
limit. As it was calculated for a probe field propagating 
freely through the condensate [9] , the depletion scales as 
the absorption which will be suppressed by choosing very 



large detuning. Therefore the cavity is essential in creat- 
ing a specific coupling of the many-atom system to the 
radiation. Because of the fast round-trips of photons, the 
cavity field experiences the collective behavior of atoms. 
This kind of coupling gives rise to the noise amplification 
mechanism analogous to the so called 'excess noise' in 
laser physics [TOj QU Q2J [13] . 

The single-mode high-Q optical cavity is driven by in- 
jecting coherent laser light through one of its mirrors. 
The laser frequency to is detuned far from the frequency 
of the atomic transition wa, i.e., |Aa| 3> 7, with the 
detuning A a = lu — lua, and 2j being the rate of sponta- 
neous emission. On the other hand, the laser is close to 
resonance with the cavity. A one-dimensional geometry 
is considered for the atomic motion. The single mode 
function of the cavity field is given by cos(fca;), k = lo/c. 
The dynamics of its annihilation operator is governed by 
the following Hciscnbcrg-Langevin-type equation of mo- 
tion in the rotating wave approximation: 

d r f - 1 

i — a(t)= -A c + / * 1 '(x,t)U(x)$(x,t)dx - in a(t) 

+ iri + i£(t). (la) 

The parameter Ac = uj — ujq is the cavity detuning, r\ 
is the strength of the cavity pumping, 2k is the photon 
escape rate, U(x) — Uocos 2 (kx) describes the frequency 
shift caused by a single atom, with U$ being the strength 
of this coupling. The noise operator £(i) describes white 
noise with the only non-vanishing correlation function 
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The equation of motion of the atomic field operator is 
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In order to distil the effect that the noise source £ in- 
filtrates into the matter wave component and induces 
depletion, the atom-atom s-wave interaction is set ar- 
tificially to zero. This is a consistent approach since 
the cavity-induced depletion is independent of the col- 
lisions. The two processes rely on distinct transitions 
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between the Bloch states of the A/2 periodic potential 
U(x). The atom-photon interaction does not couple 
Bloch states with different quasi momenta. We can thus 
consider wavefunctions being also periodic with A/2, the 
zero quasi-momcntum Bloch state representing the entire 
band. We focus on the depletion due to such 'band-to- 
band' transitions. By contrast, collisions induce transi- 
tions mostly within a Bloch band, and for usual densi- 
ties this does not influence significantly the band-to-band 
processes. The validity of the model is restricted to the 
parameter regime where the Bose gas loaded in the opti- 
cal lattice created by the cavity field is in the superfluid 
phase, far from the transition point to the Mott phase 
[HI [15] . In this regime, the small amount of depletion of 
the condensate into higher quasi-momcntum states due 
to s-wave scattering can be estimated by the usual for- 
mula well-known from Bogoliubov-theory 16J. 

The field operators are decomposed according to a(t) = 
a(t) + 5a(t) and $(x,t) = \/N(p(x,t)+8$(x,t), with a(t) 
being the coherent component of the radiation field, and 
similarly tp(x,t) stands for the condensate wavefunction 
normalized to unity, while TV is the number of atoms in 
the condensate. Substituting these into Eqs. ([!]) and ne- 
glecting those terms which are not pure c-numbers leads 
to a coupled pair of nonlinear equations similar in spirit 
to the Gross-Pitaevskii equations: 



d 



(t) = [- A c + N(U) -ws]a(t) + 



ir), 



(2a) 



h 2 A 
2m 



h\a(t)\ 2 U(x) 



<p(x,t). (2b) 



Via the average (U) — J cp* (x,t)U(x)(p(x,t)dx, the con- 
densate wavefunction ip(x) couples into the evolution of 
the field amplitude. The solution of Eqs. ^ yields the 
possible steady-states of the BEC-cavity system. 

The Heisenberg-Langevin equations (JIJ can be lin- 
earized in the fluctuations Sa(t) and 8^>(x,t) around the 
steady state of Eqs. The procedure is similar to 

that used in Ref. |17j and, for a different geometry in 
Ref. [18], however, here we will consider the effect of 
quantum noise. One obtains a coupled set of differential 
equations for the fluctuations with couplings to their her- 
mitian adjoints. For a compact notation the fluctuations 
are arranged in a vector R = [8a, 8cv, 8$>(x), 8^'(x)] for 
which the linearized equation of motion, 



d . - 
i—R = MR + i£. 
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where £ = [£, 0, 0], and the matrix M is 
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FIG. 1: (Color online) The real and imaginary parts of the 
eigenvalues of the linear stability matrix M vs. the one atom 
light shift {To- Only those eigenvalues are plotted which have 
a nonnegative real part. Every eigenvalue is degenerate at 
Uo = 0, and the lift of the degeneracy can be resolved only 
close to Uo ~ 1 for the low-lying levels, where one of the eigen- 
value in the quasi-degenerate pair is plotted by thin, dashed 
line. These are not coupled to the radiation mode and have 
vanishing imaginary part. The highest plotted eigenvalue at 
Uo = (thick, dashed line in the top panel) corresponds to 
a dominantly radiation field excitation; notice that its imagi- 
nary part is at —k. As \ Uo\ grows this mode crosses the other 
excitation energies and reaches zero at ./Wo ~ Ac- The pa- 
rameters are Ac = — 1000 ^r, k = 100 ur, r\ = IOOOojr, 
N = 1000. 



with A = -A c +N(U)-iK, Y = VN<f(x)U(x), Xf(x) = 
J ip(x)U(x)f(x)dx&ndH = -hA/(2m) + \a\ 2 U(x). The 
key point with respect to the depletion of a BEC is that 
the matter field fluctuation S^>(x) couples to both 5a and 
8a) , this latter being driven by the 'noise creation' opera- 
tor £T . L e t us also note that the matrix M is non-normal, 
i.e., it doesn't commute with its hermitian adjoint. As a 
consequence its eigenvectors are not orthogonal to each 
other, hence the analogy to excess noise in lasers and to 
the Petermann-factor arises [EJ [20j [2TJ [22] . The stability 
can be properly characterized in terms of pseudo-spectra 
by using the theory of non-normal operators [23] |2~4"] . 
However, we will restrict this work to the calculation of 
the depletion of the BEC. 

The eigenvalues u>k and the corresponding left and 
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right eigenvectors, and A k \ respectively, of M 

are calculated numerically. First the coupled Gross- 
Pitaevskii equations Q are solved on a 200 point grid 
with imaginary time propagation. From the wavcfunc- 
tions, the matrix M can be formed and then diagonalized 
numerically. The spectrum of fluctuations is exhibited in 
Fig. [I] At Uq = 0, the real part (top panel) renders the 
spectrum of a homogeneous ideal gas. In a large region 
< \Uq\ ^ 0-8 the solid lines corresponding to motional 
excitations are constant, which indicates that the con- 
densate wavefunction remains close to the homogeneous 
one. There is a critical point where the condensate pulls 
the frequency of the cavity mode into resonance with the 
pump laser frequency (N(U) = Ac), marked by that 
the real part of the eigenvalue of the dominantly cavity- 
type excitation (thick, dashed line) crosses zero. At this 
point, the steady-state wavefunction is localized in the 
dipole potential (its depth is about 100 ujr for the pa- 
rameters of the figure). Due to the finite spread of the 
wavefunction around the antinode, the resonance occurs 
for Uq slightly shifted from Ac /N. 

Coupling of the motional modes to the radiation mode 
is signified by the appearance of non-vanishing imagi- 
nary parts (Fig. [I] bottom panel). Above the critical 
point positive imaginary parts appear indicating a dy- 
namical instability of the steady-state solution. This is 
in accordance with the semiclassical theory of the me- 
chanical forces on atoms in a cavity [25] . which states 
that the cavity field heats up the motion of atoms in the 
regime of Ac — N(U) > 0. We will disregard this regime 
and consider only fluctuations around stable steady-state 
mean field solutions. Note that the clear maximum of 
the imaginary parts reflects the semiclassical result that 
the cavity cooling mechanism is most efficient at half a 
linewidth below resonance [25] . 

The left- and right eigenvectors of M can be used to ex- 
pand the fluctuation vector R in terms of "quasi-normal" 
modes: R = J^fcPfcE • By use of the orthogonality of 
the left and right eigenvectors, (l^ k \r^) — Sk,i, where 
(a,b) is the Euclidean scalar product, the quasi-normal 
mode amplitude operators are obtained as p k — (l^jR)- 
They evolve independently 



expressed in terms of the quasi-normal modes as 
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where the projected noise is Q k = (l (k \i). With the inte- 
gration of Eq. ([5]) the time dependence of the fluctuation 
operators read 



(t) = e-^pkiO) + / e-^-^Q^Odi'. (6) 



Starting from the definition of the depletion, it can be 



SN(t)= j ^s¥(x,t)S^(x,t)^dx 

= X>*(*)M*)> / rf\x)r<i\x)dx, (7) 
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with 



l_ e -i(^,)t , 

• f \ \ 2 

l{U! k + Uli) 



(8) 



By assuming that we start with a pure BEC at t = 0, the 
first term on the right-hand side of Eq. Q is zero. The 
depletion then reads 

5N{t) =2 K J2 ., - — A ] 4° 



i(to k + lui) 



x / r[ k \x)r ( 3 l \x)dx. (9) 
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Contrary to the case of a conservative system, where 
quantum depletion is the sum of the individual contri- 
butions of the noncondensed one-particle state occupa- 
tion numbers, here the depletion is the double sum of 
contributions from pairs of quasi-normal modes (k and 

0- 

The matrix M has an important symmetry property, 
namely rMX = — M*, originating from the fact that 

R^ = TR, with r being the permutation matrix that ex- 
changes the first row with the second and the third row 
with the fourth one. This property ensures that if uj k 
is an eigenvalue, — lu^ will also be an eigenvalue. There- 
fore the denominator in Eq. (JoJ) may become small for a 
pair of eigenvalues (u> k and u>i) if He{uj k } — — Re{w;}, i.e. 
when the quasi-normal modes k and I are connected with 
the transformation T. Moreover, the smaller the imagi- 
nary parts (their dampings), the smaller the denomina- 
tor. However, an eigenvalue with a small imaginary part 
is accompanied by an eigenvector with small "photon" 
components (diminishing l\ and I2 terms). Therefore not 
all of the eigenvalues with small imaginary parts play a 
significiant role and the evaluation of the terms remain a 
numerical task. Nevertheless, the contribution to the de- 
pletion from such pairs of quasi-normal modes dominates 
the double sum. 

The total steady-state depletion in the t — > 00 limit 
[the exponential in Eq. Q is dropped] is shown in Figure 
[2] (semilog scale). The thick solid line corresponds to the 
parameter set of the spectrum in Fig. [T] The depletion 
is closely constant for a wide range of Uq up to the reso- 
nance regime defined by | Ac — N(U)\ < k. At resonance 
the depletion diverges, in accordance with the fact that 
the dynamical equilibrium of the system becomes unsta- 
ble. It is interesting that there is a dip slightly below the 
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tical cavity is exposed to excess quantum noise which 
indicates that the long-time evolution drives the system 
into exotic quantum phases unexplored so far. 

We acknowledge funding from the National Scientific 
Fund of Hungary (NF68736, T046129, and T049234). 



0.2 0.4 0.6 0.8 1 

\Uq\ in units of lj r 

FIG. 2: (Color online) The depletion SN of the condensate as 
a function of the coupling constant Uo for various detunings. 
The parameters are k = 100 ujr, r\ = —Ac, and N = 1000. 
In the plotted range of Uo, for smaller detunings (solid and 
dotted lines) the BEC-cavity system can be resonant with 
the pump, which leads to divergences. Well below resonance, 
iV|Z7o| <C Ac, the depletion depends weakly on Uo (see the 
plateau of the solid and the dashed lines) and on N. In this 
regime the depletion scales with Ac/cor. 



resonance. We attribute it to the significant change in 
the shape of the condensate wavefunction, from homoge- 
neous to a strongly localized one, which yields a variation 
of the overlap factors in the matrix M. It is important to 
note that the finite, large value of depletion in the limit 
of Uo — *■ was calculated by taking first the t — > oo limit 
of the time dependent result Eq. ([£]). However, during a 
finite measurement time the system does not necessarily 
reach the steady state since the relaxation time diverges 
as Uq 2 for small Uo. 

The actual steady-state value of the depletion is deter- 
mined by the ratio of the cavity mode frequency Ac and 
the recoil frequency uir, i.e., the relation of the energy 
scales describing optical excitations and atomic motion, 
respectively. This is illustrated by the other two curves 
in this figure. Dashed line corresponds to a detuning Ac 
increased by an order of magnitude, then the depletion 
is also scaled up by a factor of 10. Obviously, the res- 
onance regime is pushed out of the plotted range and 
only the initial plateau can be seen. Oppositely, when 
decreasing Ac by a factor of 10, the initial depletion de- 
creases (dotted line). However, for such a small detuning 
Ac ~ K the system is in the resonance regime already 
for small Uq values, the wide plateau is missing and the 
divergence of the depletion is exhibited. As long as the 
atom-cavity system is far from resonance (on the wide 
plateau), the depletion SN is independent of the atom 
number N. Similarly, the photon number \a\ 2 is quite 
irrelevant to the amount of depletion, it influences only 
the shape of the ground state in the optical potential and 
thereby some overlap factors in the matrix M. 

In conclusion, we proved that a BEC in a high-Q op- 
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